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Abstract 

In all spatial dimensions d, we study the static and dynamical properties of a general- 
ized Smoluchowski equation which describes the evolution of a gas obeying a logotropic 
equation of state, p = Alnp. A logotrope can be viewed as a limiting form of polytrope 
(p = Kp" 1 , 7 = 1 + 1/n), with index 7 = or n = -1. In the language of generalized 
thermodynamics, it corresponds to a Tsallis distribution with index q = 0. We solve 
the dynamical logotropic Smoluchowski equation in the presence of a fixed external force 
deriving from a quadratic potential, and for a gas of particles subjected to their mutual 
gravitational force. In the latter case, the collapse dynamics is studied for any negative 
index n, and the density scaling function is found to decay as r~ a , with a = for 
n < -|, and a = for -| < n< 0. 



Key words: Nonlinear meanfield Fokker-Planck equations, generalized thermodynam- 
ics, polytropic equation of state, self-gravitating systems, chemotaxis. 

1 Introduction 

Recently, several researchers have studied generalized forms of Fokker-Planck equations QjJ E] . 
These equations arise when the coefficients of diffusion, mobility and friction in the usual 
Fokker-Planck equations (Kramers, Smoluchowski,...) explicitly depend on the concentration 
of particles. This can take into account phenomenologically "hidden constraints" that are not 
directly accessible to the observer jSj. These nonlinear Fokker-Planck equations are associated 
with generalized free energy functionals and non-standard distributions. For example, Tsallis 
g-distributions constitute an important class of non-Boltzmannian distributions 0]. Plastino 
& Plastino [5] have shown that they could be obtained as stationary solutions of a nonlinear 
Fokker-Planck equation, with potential applications to the context of porous media. Recently, 
Kaniadakis 6], Frank [7] and Chavanis [H] have shown that it was possible to construct more 
general Fokker-Planck equations leading to other forms of non-standard distributions. For 
example, Chavanis jHl El ED] has introduced a generalized Smoluchowski equation involving 
an arbitrary barotropic equation of state p = p(p). This can be derived from a generalized 
Kramers equation in a strong friction limit jHl • The usual Smoluchowski equation leading 
to the Boltzmann distribution corresponds to an isothermal equation of state p = pT while the 
Tsallis distributions are associated with a polytropic equation of state p = Kp 1 where 7 = 1 + ^ 
plays the role of the Tsallis g-parameter. 
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In this paper, we consider the case of a logotropic equation of state p = A In p where A 
is a positive constant. This equation of state sometimes appears in astrophysics 12] > but we 
consider it here at a more general level. In Sec. 121 we show that a logotrope can be viewed 
as a limiting form of a polytrope with 7 = (or as a Tsallis distribution with q — 0). In 
Sec. El we solve the logotropic Smoluchowski equation with a fixed quadratic potential. This 
is a particular case of the general analytical solution found by Plastino & Plastino (3] and 
Tsallis & Buckman ^3] fo r the nonlinear Fokker-Planck equation (associated with polytropes). 
In Sec. |H we consider a coupling with the gravitational potential [HI [T5] and study the 
logotropic Smoluchowski-Poisson (LSP) system. This is a particular case of the nonlinear 
(polytropic) Smoluchowski-Poisson system studied in jTHj, but with a negative index n = — 1 . 
We characterize the stationary solutions of the LSP system, determine their stability and find 
self-similar collapse solutions when no equilibrium state exists (in a space of dimension d > 2). 
We find that the collapse of polytropes with negative index is peculiar: (i) for n < —d/2, the 
scaling exponent is a = 2n/(n — 1) like for positive index (ii) for — d/2 < n < 0, the scaling 
exponent is a = 2d/(d + 2) corresponding to the collapse of a cold (pressureless) system. In 
particular, the scaling exponent of logotropes is a = 2d/(d + 2) (d > 2). 



2 Nonlinear mean-field Fokker-Planck equations 



2.1 The generalized Smoluchowski equation 

In a recent series of papers (13 US1 U3 UH UH EH EU 1221 122] , we have studied a general class 
of nonlinear mean-field Fokker-Planck equations of the form 



(1) 



dt 



- (Vp + pV$) 



(2) $M) = J u(\r-r'\)p(r',t)dr', 

introduced in Chavanis [HI HH CHI 121] • This type of equations arises in different domains of 
physics, chemistry, astrophysics and biology j2H]- The drift-diffusion equation (JH) can be viewed 
as a generalized Smoluchowski equation where p = p(p) is a barotropic equation of state and 
$(r, t) is a potential. This potential can be imposed by an external medium (in which case 
$ = $ ex t(r) is assumed given) or generated by the particles themselves through the mean-field 
equation (j21) where u(\r — r'|) is a binary potential of interaction. The generalized Smoluchowski 
equation can be derived from a generalized Kramers equation in a strong friction limit 
£ — > +00 by performing a Chapman- Enskog expansion in terms of the small parameter l/£ 
|TT] or by using a method of moments [TOj- It can also be directly obtained from the damped 
barotropic Euler equations by neglecting the inertial term in the equation for the velocity [S]. 
Thus, the generalized Smoluchowski equation can be interpreted as an overdamped limit of 
kinetic or hydrodynamic equations in which a friction force (real or effective) is present [2*2*] . 
The generalized mean-field Smoluchowski system CQ)-© decreases the Lyapunov functional 

(3) F[p] = Jpf ^dp' dv+ l -Jp<$> dr, 

which can be interpreted as a generalized free energy 0. Indeed, F = - J -^(Vp+pV$) 2 rfr < 0. 
The stationary solutions of Eq. ((TJ) satisfy the condition 

(4) Vp + pV$ = 0, 
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which can be viewed as a condition of hydrostatic equilibrium. Using Eq. (@J) and the equation 
of state p = p(p), we find that the equilibrium density field is of the form p = p[$(r)]. When 
substituted in Eq. (J2J), this yields an integro-differential equation for the potential $(r). The 
steady states of Eqs. (HI)-© are extrema of the free energy (j3J) at fixed mass M = f p dr. They 
satisfy the first order variations d~F — a5M = which directly return Eq. (@J). Furthermore, the 
steady states are linearly dynamically stable if, and only if, they are minima of the free energy 
(JSJ) at fixed mass M = J pdr 0. 

When $ is interpreted as the gravitational potential determined by the Poisson equation, 
the stationary solutions of the generalized Smoluchowski-Poisson system (GSP) are the same 
as the stationary solutions of the Euler- Poisson system 

(5) | + V( P u) = 0, 



(6) -^ + (u-V)u = — Vp + V$, 

at p 



(7) A$ = S d Gp, 



describing barotropic stars [2d] . However, the dynamics of these two systems is different 
as the Euler-Poisson system describes a fluid without friction £ = while the generalized 
Smoluchowski-Poisson system describes a Brownian gas in an overdamped limit £ — > +oo [2~2"] . 
With these analogies and differences in mind, we can use many results obtained in astrophysics 
[26 when studying nonlinear mean-field Fokker-Planck equations of the form We think 

therefore that it is enlightening to use a similar vocabulary and similar notations. 

We note, finally, that drift-diffusion equations of the form also arise in biology, in 

the context of chemotaxis fHlEZI- They can be written in the form 

(8) ^ = V • [x (Vp - pVc)] , 



dc 

(9) — = -k(c)c + f(c)p + D c Ac. 

They describe the collective motion of cells (usually bacteria or amoebae) that diffuse and that 
are attracted by a chemical substance that they emit. Here, p(r, t) denotes the cell density and 
c(r, t) is the concentration of chemo-attractant which induces the drift force. This model has 
been introduced by Keller & Segel j2Hj- In its most studied version [22], p(p) = (-D*/x)p is a 
linear function of the density which is similar to an isothermal equation of state. However, the 
general Keller-Segel model [2H] allows the possibility that the diffusion be anomalous so that 
the function p = p(p) can be nonlinear. Note also that in the limit of large diffusivity of the 
chemical D c — > +oo (and for k(c) = k and f(c) = f), the time derivative dc/dt can be neglected 
[22] so that Eq. (j2J) reduces to a Poisson equation (like in gravity) with a screening term —kc. 
Therefore, the general study of equations of the form Q-Q connects several active areas of 
research in theoretical physics: nonlinear Fokker-Planck equations, self-gravitating systems, 
chemotaxis and non-extensive thermostatistics. Here, we shall focus on the case where the 
equation of state is that of a logotrope p = A In p and we shall make the connection with the 
previously studied polytropes p = Kp 1 . 
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2.2 The polytropic Smoluchowski equation 

For the polytropic equation of state p = Kp 1 (with 7 = 1 + 1/n), we obtain the polytropic 
Smoluchowski equation 



(10) 



dp 
di 



V 



{KVp 1 + pV$) 



The free energy (j3J) corresponding to a polytropic equation of state can be conveniently written 



'IF 



F W = ^ZT[ |(p 7 -p)rfr + l /p$dr, 



and the stationary solutions of Eq. (JTUJl are given by 
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7-1 

A — $ 



This polytropic distribution reproduces the statistics introduced by Tsallis in his generalized 
thermodynamics j3J. When $ is a fixed external potential, the polytropic Smoluchowski equa- 
tion (flop is equivalent to the nonlinear Fokker-Planck equation introduced by Plastino & Plas- 
tino [5 . When $ = 0, it returns the equation of porous media. We note that the free energy (|llj) 
can be written F = E — T e gS with E = | f p$ c/r, T e g = K and S = — -^ f (p 7 — p)dr. In the 
language of generalized thermodynamics, this can be viewed as a Tsallis free energy with index 7 
where K plays the role of a generalized temperature [Sj. For 7 = 1, i.e. n = +00, we recover the 
isothermal equation of state p = pT and the Boltzmann free energy F = | J p$ dr+T f p In p dr. 
As we shall see, the case 7 = 0, i.e. n = —1, is also special and corresponds to what have been 
called logotropes in astrophysics [12] . 

2.3 The logotropic Smoluchowski equation 

Let us consider the logotropic equation of state 

(13) p = A\np, 

where A is a constant. This equation of state has been introduced in astrophysics to account 
for certain properties of molecular clouds that could not be understood in terms of isother- 
mal distributions ^2]- However, this equation of state can have application in more general 
situations, beyond the realm of astrophysics, so we shall consider it here on a general footing. 
Inserting Eq. (JT3J) in Eq. (JTJ), we get the logotropic Smoluchowski equation 



dp 
di 



(AV In p + pV$) 



(14) - = V • 
The free energy (J3J) of a logotrope takes the form 

(15) F{p] = -A J \npdr+^ J p$dr. 



It can be written F = E — T P frS with T P 



A and 



(16) 



S[f) 



hap dr. 
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In a previous paper [3U], this functional has been called the log-entropy. The stationary states 
of the logotropic Smoluchowski equation are given by 



< 17 > "=ATi¥' 

Note that when the potential is quadratic, i.e. $ = r 2 /2, the equilibrium distribution is the 
Lorentzian as noted in [3U] • Thus, the log-entropy (fTUJl can be seen as the functional associated 
with the Lorentzian distribution. 

2.4 The polytropic index 7 = 

Returning to the polytropic Smoluchowski equation (fTUj). we note that the index 7 = is 
special because the polytropic Smoluchowski equation (fTUj) and the Tsallis free energy (fTTj) 
become trivial for this index except ii K ^ +00 [3U]. This type of indetermination is often 
characteristic of a logarithmic behavior as can be seen with a simple change of variables. If we 
set K — A/7, the polytropic Smoluchowski equation (fTUJ) can be rewritten 



Us) | = v. 



Under this form, this equation describes polytropes for 7 7^ and logotropes for 7 = 0. 
Therefore, a logotrope can be viewed as a polytrope with index 7 = 0, i.e. n = —1. The 
stationary solution of equation (fTSj) is given by 



(19) p 



7-1 

A — — — <3> 

A 



7-1 



and this distribution passes to the limit for 7 — > 0. Indeed, for 7 7^ we can put Eq. (fT9"j) in the 
form of Eq. (fT2~j) and for 7 = 0, we recover Eq. (fTTj). The free energy associated with Eq. (fTHj) 
can be written 

(20) F[p] = —A- J(p'r-p)dr+ 1 -J p$ dr. 

For 7 7^ 0, we can put Eq. (|2Uj) in the form of Eq. (fTTj) and for 7 = 0, we recover Eq. (fTB]) up 
to additive terms that do not depend on the density. Since only the variations of F matter, 
these terms (that can be infinite!) do not play any role (3H|. In the same spirit, the general 
polytropic equation of state can be written p = ^p 7 with the convention that lim 7 _ > o — = 
lnp + Cst. Therefore, in the language of generalized thermodynamics, logotropes correspond to 
Tsallis statistics with index 7 = 0. 

2.5 Other formulations 

We can write the generalized Smoluchowski equation in the form 

(21) ^ = V-[D(pC"(p)Vp + /3pV$)], 

where C(p) is a convex function, (3 is a constant which can be interpreted as a generalized 
inverse temperature and D is a positive constant which can be interpreted as a generalized 
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diffusion coefficient. The mobility p, = l/£ satisfies a form of Einstein relation p, = Df3. From 
the convex function C we define a generalized entropic functional 



(22) S[p] = - J C(p)dv, 

and a generalized free energy J = S — j3E. This functional increases monotonically with time 
(J = f — (pC"(p)Vp + /3pV$) 2 <ir > 0). This can be viewed as a form of if-theorem in a 



canonical description where the inverse temperature (3 is fixed jSJ. The steady solutions of 
Eq. (|21|) are determined by the integro-differential equation 

(23) C'(p) = - a, 

where $ is related to p by Eq. (J2J) in the general case. They are critical points of J at fixed mass 
M; indeed, the first order variations satisfy 5J—a5M = which returns Eq. (J23|) . Furthermore, 
the linearly dynamically stable stationary solutions of Eq. ()21|) are maxima of J at fixed mass 
M. 

For the Tsallis entropic functional written in the form 

(24) 5W= __L_y (P 7_ p)drj 

the generalized Smoluchowski equation (J2"T|) reads 

(25) !| = V ' [Dip'-'Vp + PpV®)}. 

Comparing with Eq. (jl8j) . we find that (3 = I /A. This is consistent with the previous formalism 
where A = T e ff is interpreted as an effective temperature. The stationary solution of Eq. (|23j) 
can be written 

(26) p= [X-P(j- 1)$]^ . 

We note that the parameter (3 has not the dimension of an inverse temperature, although this 
is the quantity which naturally enters in the free energy functional J and in the variational 
principle SS — [35 E — a5M = determining the stationary solution (jSHJ). As discussed in |3*T] . 
there are different notions of temperature in the case of polytropic distributions. We can write 
the distribution (J26|) in the alternative form 

(27) p = g[l_ fo ( 7 _l)$]^T ; 

where Q = X~ and b = (3/X now has the dimension of an inverse temperature. 



2.6 The logotropic Kramers and Landau equations 

The generalized Fokker-Planck equations (JIJ and ()21j) are written in position space. We can 
also introduce generalized Fokker-Planck equations in velocity space. When the friction term 
is linear in v, they are referred to as generalized Kramers equations. Some exemples are 
given in [B] for different forms of entropic functionals S = — J C(f) drdv. For the log-entropy 
C(f) = —Inf considered in this paper, we obtain the logotropic Kramers equation 



(28) 



df = d_ 
dt dv 
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The generalized Kramers equation is associated with a canonical description where the tem- 
perature (3 is fixed and the free energy J = S — [3E increases [HJ 121] • Alternatively, we can 
consider the generalized Landau equation that is associated with a microcanonical description 
where the energy E = J f^-dvdv is conserved and the generalized entropy S increases 
For the log-entropy C(f) = — In/, we obtain the logotropic Landau equation 



(29) 



?L = a— 

dt dv v 



fl d\nf f d\nf 



dv 1 



dv' 1 



dv', 



where / = /(v, t) and /' = /(v', t). The steady solution of the logotropic Kramers and Landau 
equations is the Lorentzian 



(30) 



fe 



1 



A + /3f 

Note that if the velocity is not bounded, this distribution is not normalizable in d > 2. 



3 The logotropic Smoluchowski equation with a fixed 
quadratic potential 

We shall here provide an analytical time- dependant solution of the logotropic Smoluchowski 
equation in d — 1 for a fixed quadratic potential $ = x 2 /2 + Ax. This solution is a straight- 
forward extension of the solution found by Plastino & Plastino [H] and Tsallis & Bukman ^H] 
for the polytropic Smoluchowski equation with 7 ^ 0. We cannot directly make 7 = in their 
solution for the reason discussed in Sec. 12.41 However, if we write the polytropic Smoluchowski 
equation in the form 



(31) 



dp d 
dt dx 



ox 



we can follow exactly the reasoning of EI] and finally let 7 — > since this limit is now well- 
behaved. Note that Eq. (j31|) with A = coincides with the polytropic and logotropic Kramers 
equation if x is replaced by v jHJ. 

The stationary solution of Eq. (}3l"| can be written 



(32) 



P = Q 



-6(7-1) (x + 



1 

7-1 



where b = fi/Q 1 ^ 1 and Q is determined by the normalization condition (see below). For 7 > 1, 
the density goes to zero at finite values of x — xi, x 2 (say) and it is implicitly understood that 
p = for x < X\ and x > x<i- For 7 < 1, the density distribution is defined for all x and it 
decreases algebraically as x -2 ^ 1 " 7 - 1 . Following the original idea of Plastino & Plastino |SJ, we 
look for time-dependant solutions of the polytropic Smoluchowski equation in the form 



(33) 



p(x,t) = Q(t) 



2*(*)(7 



l)(x-m(t)f 



7-1 



where Q(t), b(t) and m(t) are functions of time, 
requires that 7 > — 1 and leads to the relation 



The normalization condition f_ °° pdx = 1 



(34) 



b(t) = 2Q(t)V 



7 



where we have defined 

(35) // = / [1 - (7 - l)y 2 ] i- 1 dy. 

J — oo 

Substituting the Ansatz ()33|) in Eq. (|3ip. we find after straightforward calculations that 

(36) m = —(3m — A, 



(37) 



Q = -2fi 2 Q-* +2 + (3Q. 



Starting from a Dirac distribution p(x, 0) = <5(x — xo), these equations are readily integrated 
and we obtain 



(38) 



Q(t) 



P 



2p? 1 - e -^ +1 )* 



7+1 



(39) 



(40) 



b{t) = 2p 2 Q(t) 2 



For 7/O, this returns the results obtained in Ej. However, we can now pass to the 
logotropic limit 7^0 and get 



(41) 



p(x,t) 



Q(t) 



l + \b{t){x-m{t)f 



(42) 



Q(t) 



(3 1 
2tt 2 1 -e-/ 3 *' 



(43) 



m(0 = zoe- /9t + ~(e-' 9t -l), 



(44) 6(t) = 27r 2 Q(t) 2 . 

In particular for /3 = 0, we obtain 

1 



(45) 
Hence 
(46) 



Q(t) = ^7' m (*) = x o - At, b(t) 



2n 2 t 



2vr 2 t 2 ' 



p(x,t) 



27T 2 t I + (x-xo) 
1 ^ 4ir 2 t 2 
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4 The logotropic Smoluchowski-Poisson system 



We shall consider here the coupling between the logotropic Smoluchowski equation and gravity. 
Thus, we consider the logotropic Smoluchowski-Poisson (LSP) system 



(47) 



dp 
dt 



(AV lnp + pV$) 



(48) A$ = S d Gp, 

in d-dimensions 1 . This can be viewed as a particular case of the polytropic Smoluchowski- 
Poisson system studied in corresponding to the index 7 = (n — — 1). This case must 
be treated separately (i) because the study of [T3] is restricted to n > and (ii) because the 
index 7 = needs a special discussion as we have seen previously. As indicated in Sec. I2.1| the 
steady solutions of the LSP system are the same as the steady solutions of the Euler-Poisson 
system with an equation of state p = A In p describing logotropic stars and logotropic clusters 
in astrophysics [TJ] . Therefore, the results of Sees. I4.1M.7I for the steady states also apply 
to logotropic stars and logotropic clusters, while the results of Sec. 14.81 for the dynamics only 
apply to the LSP system. 



4.1 The logotropic Lane-Emden equation 

For a spherically symmetric distribution, the equation of hydrostatic equilibrium (0J) combined 
with the Gauss theorem d^/dr = GM(r)/r d ~ l can be written 

(49) * = -„ GM M 



dr r d 1 ' 

where M(r) = J p(r' \t)Sdr' d ~ 1 dr' is the mass within the sphere of radius r. The foregoing 
equation can be put in the form 



(50) ^ d fr d 1 dp \ g ^ 

nr.d-1 dr \ p dr 



For the logotropic equation of state p = A In p, we define 

(si) P = Po /e, £=(^) 1/2 ^ 

where po is the central density and £ a scaled distance. Substituting these relations in Eq. ([HU|. 
we obtain what might be called the logotropic Lane-Emden equation 

(52) (e"- 1 -) - - 

with 6(0) = 1 and 0'(O) = 0. 



1 It should be made clear that, in the present work, we consider a self-gravitating Langevin gas with a 
prescribed polytropic or logotropic equation of state. We furthermore consider a strong friction limit in which 
the dynamics is governed by the (generalized) Smoluchowski equation. This context is very different from the 
works of Taruya & Sakagami (141 132) who consider the Hamiltonian TV-stars system described by the orbit- 
averaged-Fokker-Planck equation and show, by means of direct numerical simulations, that the transient phases 
of the dynamics (quasi-equilibrium states) can be fitted by a sequence of polytropic distributions (Tsallis) with 
a time dependant index n(t) or q(t). 
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4.2 The singular logotropic sphere 



We note that for d > 1, the function 



(53) 9 - = ^ 

is solution of the logotropic Lane-Emden equation. This is, however, not a regular solution since 
it does not satisfy the boundary conditions at £ = 0. In terms of the density, this corresponds 
to a singular sphere 



.4 



nl/2 



(54) [S d (d-1)G_ 

whose density diverges at the origin. It is interesting to compare this singular logotropic sphere 
whose density profile decreases as r -1 for d > 1 to the singular isothermal sphere p s = [2(d — 
2) / (SdGj3)]r~ 2 whose density profile decreases as r~ 2 for d > 2 [T2j. We note that the total mass 
of the system diverges for r — > +oo. However, if we consider box-confined configurations (in a 
sphere of radius R), the density distribution of the singular logotropic sphere can be written 



(55) 



Ps 



M(d- V 
S d R d 



-i 



It should be finally recalled that such singular solutions are not limited to isothermal and 
logotropic distributions: polytropic spheres with index n > n 3 = d/(d — 2) (for d > 2) and 
n < — 1 (for d > 1) also admit singular solutions scaling like r~ a with a = 2n/(n — 1) [To] . 



4.3 Logotropic profile in d = 1 

In d = 1, the logotropic Lane-Emden equation becomes 

(56) dv = e> 

with 9(0) = 1 and #'(0) = 0. This is similar to the equation of motion of a particle in a potential 
V(6) = — ln.6 where 9 plays the role of position and £ the role of time. Using the boundary 
conditions at £ = 0, the first integral is 

(57) Kir^ ' 
leading to 

Performing the change of variable ln0 = y 2 , we obtain the equivalent expression 

(59) V2 / e y2 dy = f . 

Jo 

For £ — ► +oo, we get 



(60) 6 ~ Inf. 

Therefore, the density profile decreases like l/(rVmr). The mass within a domain of size r 
behaves as M(r) ~ (lnr) 1 / 2 for r — > +oo so that the total mass in infinite. The density profile 
of a one-dimensional logotrope is represented in Fig. ^ 
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Figure 1: Density profile of a logotrope in d — 1. 
4.4 Asymptotic behaviors 

We shall determine the asymptotic behaviors of the solutions of the logotropic Lane-Emden 
equation For £ — > 0, using a Taylor expansion, we obtain 



(61) 



f = 1 + 7rA 



2eT 8ci(rf + 2) 



£ 4 + - ce-o). 



To investigate the limit £ — > +oo, we first perform the change of variables t = ln£ and 9 = £z. 
Equation ([52)1 then becomes 



(62) 



<i 2 ^ dz 1 



This is similar to the damped motion of a particle in a potential V(z) = -lnz + (d - l)z 2 /2. 
For d > 1, the potential has a minimum at z c = \j\Jd — 1. Thus, for t — > +oo the "particle" 
will reach this minimum, i.e. 2 — > z c . Returning to original variables, this implies that # — > 
# s = £/ 'y/d— 1 for £ — > +oo. To get the next order correction, we set z = z c + z' where z' <C z c . 
Keeping only terms that are linear in z f , we obtain 



(63) 



d 2 z' dz' 



v d 2 -2d + 2 
d- 1 



z' = 0. 



Noting that the discriminant A(d) = d 2 — 4(d 2 — 2d + 2) 2 / (d — l) 2 of the associated quadratic 
equation is always negative for d > 1, we find that 



(64) 



= e- d '/ 2 cos( + . 



Returning to original variables, we obtain the asymptotic behavior 
(65) 0~0 s (£)(l 



C /A, , r 
cos I In 4 + o 



£d/2 



(£^+oo). 
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Figure 2: Density profile of a logotrope in d = 2 and d = 3 (the oscillations are hardly visible 
but they exist). 

As in the case of isothermal and polytropic distributions |T71 ITS] , the density profile of a 
logotrope presents damped oscillations around the singular solution. The density profiles of a 
logotrope in d = 2 and d = 3 are represented in Fig. |21 For these dimensions A (2) = —12 and 
A(3) = -16. 



4.5 The Milne variables 

As in the case of polytropes [T5|, it is convenient to introduce the variables 



(66) 



u 



99' 



£_ff_ 
9 ' 



which are the appropriate forms of the Milne variables [SB] in the present context. Using the 
logotropic Lane-Emden equation (|52jl we easily derive 



(67) 



1 du 1 . 
ud£ £ 



(68) 



1 dv 1 , n . s 
-- = -2-d + u-«, 
v d£ £ 



so that the function v(u) satisfies the first order differential equation 

u — v — d + 2 



(69) 



u dv 
v du 



u + v — d 



As for isothermal and polytropic spheres, the reduction of a second order differential equation 
(Jo"2|) to a first order differential equation is a consequence of the homology theorem The 
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Figure 3: The solutions of the logotropic Lane-Emden equation in the Milne plane (u,v) for 
d — 1, d — 2 and d = 3. 



solution curve in the Milne plane is parameterized by £ going from to +oo. It is represented in 
Fig. Eland forms a spiral for d > 1. The points of horizontal tangent correspond to 2—d+u—v = 
and the points of vertical tangent to d — v — u — 0. They coincide for u s — d — 1 and t> s = 1 
which corresponds to the singular logotropic sphere ()53|) . On the other hand, the curve starts 
at (d, 0) for £ = with the slope (dv/du)o = —(d + 2)/d and rolls up to the singular logotropic 
sphere (when d > 1) for £ — > +oo. For d — 1, the solution curve tends to (0, 1) for £ — > +oo. 
It does not make a spiral but it presents a maximum at some point. 



4.6 The series of equilibria 

Since the density profile of logotropic distributions in d > 1 decreases as r _1 at large distances, 
the total mass is infinite (the total mass is also infinite in d = 1). Therefore, we shall enclose 
these distributions in a spherical box of radius R as in the case of isothermal and polytropic 
spheres |17| ITB] . For box-confined logotropes, the solution of Eq. (|52|) is terminated by the box 
at a normalized radius 



SdGpl 



1/2 



(70) a ={^) K 

The total mass of the configuration can be written 

It is therefore natural to introduce the dimensionless parameter 

/ r \l/2 i 

(72) Tj = M 



S d A J R d ~ 
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Then, the series of equilibria is defined by 



(73) 



n = e'(a). 



This relation can be interpreted in different manners. For example, if we fix A and R, the 
parameter a is proportional to the central density po and the parameter rj to the mass M. 
Therefore, Eq. ([75)1 can be interpreted as the mass-central density relation. Alternatively, if we 
fix M and R, then the parameter rj is related to A = T e ff which is interpreted as a "generalized 
temperature" . Thus Eq. (J73j) gives the relation between the central density and the generalized 
temperature. In terms of the Milne variables, the relation (|T3*|) can be rewritten 



where, by definition, uq = u(a) and vq = v(a) are the values of the Milne variables at the box 
radius. The series of equilibria ([73~)l is plotted in Fig. |U and it presents damped oscillations for 
d > 1 (for d = 1 the curve is monotonic). Similar oscillations are encountered for isothermal 
spheres in Newtonian gravity [31] and general relativity [35] ■ For a — > +oo, we recover the 
singular sphere with rj s = 1/yd — 1. The control parameter rj(a) is extremum for a c such that 
drj/da = 0. Using Eqs. (jfiT)) - ([()%]) . we find that this condition is equivalent to 

(75) u = d — 1 = u s . 

The number of extrema can be obtained by a simple graphical construction in the Milne plane. 
They are determined by the intersection between the (u, v) curve and the line u = u s . Since the 
straight line u = u s passes by the center of the spiral, there is an infinite number of extrema. 
The critical value of 77 corresponding to the first maximum is 7] c = 0.715657... in d = 3 and 
1] C = 1.028728... in d — 2. Clearly, there is no steady state for rj > r] c . Similar results are found 
for isothermal systems j3UE[- Note, however, that the series of equilibria i](a) for logotropes 
already presents oscillations in d = 2 unlike isothermal spheres for which oscillations appear 
for d > 2. 

4.7 Stability analysis 

We shall now study the dynamical stability of logotropic spheres by determining whether they 
are minima of the generalized free energy functional (fTo^) at fixed mass. As shown in [H], a sta- 
tionary solution of the generalized Smoluchowski-Poisson system ©-(0) is linearly dynamically 
stable if, and only if, it is a minimum of the functional F defined by Eq. (j3]). In that Brownian 
context, F is interpreted as a Lyapunov functional or as a (generalized) free energy, and its 
minimization corresponds to an effective thermodynamical stability criterion in the canonical 
ensemble. On the other hand, a stationary solution of the barotropic Euler-Poisson system 
©-© is nonlinearly dynamically stable if it is a minimum of F. In that Euler context, F is 
related to the energy functional of the barotropic gas W (see [I5II31II35] for details). Therefore, 
our stability analysis has applications for these two different systems: a Brownian gas described 
by nonlinear Fokker-Planck equations and a "normal" gas described by the Euler equations. 
Since the following stability analysis is very similar to that developed in [T5] for polytropic 
distributions, we shall only give the main lines of the analysis and refer to ^5] for more details. 
The second order variations of F are given by 



(74) 




(76) 
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0.4 0.6 
r/R 



Figure 5: Perturbation profile corresponding to the first mode of instability at r\ c in d = 2 and 
d = 3. 



15 



Restricting ourselves to spherically symmetric perturbations, introducing the quantity q(r) 
through the defining relation 

^ ^ ^ Sd,r d ~ l dr' 

and using integrations by parts, we can put the second order variations of F in the quadratic 
form 



( 78 ) PF = —I ^[-a^r^ + ^Qdr. 



1 r R d f A d G 

2 Jo ^ dr \SdP 2 r dr ~ 1 dr r 

We are thus led to considering the eigenvalue equation 



(79) — ( A - | G ^ q (r) - Xq (r) 

dr \SdP 2 r d ~ l dr r d ~ l J 

and determine when, in the series of equilibria, the eigenvalues pass from positive (stable) to 
negative (unstable) values. To that purpose, it is sufficient to determine the point of marginal 
stability corresponding to A = 0. Introducing the dimensionless quantities of Sec. 14.11 we have 
to solve 

d f 6 2 dF" 

dt Ve^^y 1 e 



(8°) ^ = ^(^5=T^ +75=r = °. 



with the boundary conditions F(0) = and F(a) = implied by the conservation of mass. 
Noting that 



£(£ d -V) = 29', c(tp)=2(d-l)&, 



we find that the solution of Eq. (fKUj) satisfying F(0) = is 
(82) F = 



^ - (d-i^e' 





where a is an arbitrary constant (due to the linearity of the eigenvalue equation). The critical 
values of a where an eigenvalue becomes zero are determined by the condition F(a) = 0. 
Using Eq. (|82|) and introducing the Milne variables we find that this is equivalent to 

uq = d — 1 = u s . Comparing with Eq. (j75j) . we find that a new eigenvalue becomes zero each 
time that rj is extremum. As a result, the series of equilibria becomes unstable at the point of 
maximum 7], previously denoted r) c . This result is consistent with the turning point argument 
of Poincare (see ^5] for details). The perturbation triggering the instability at the marginal 
point is given by 

po S d £ d 1 d£ 

In terms of the Milne variables, this can be rewritten 
(84) ^ = f («.-.)■ 

Hence, the number of oscillations of the perturbation profile 5p(£) can be easily determined by 
graphical constructions. The values of £ where 5p(£) vanishes are determined by the intersection 
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between the (u, v) curve and the line v = v s . Of course, we have to consider only values of 
£ < a. At the point of instability t] c , the perturbation has only one node (see Fig. EJ). 

We can also study the linear dynamical stability of a stationary solution of the LSP system 
and show the equivalence with the stability criterion based on the minimization of the free 
energy. Linearizing the generalized SP system around a stationary solution and writing the 
perturbation in the form Sp ~ e A< we finally arrive at the eigenvalue equation jH] : 

(85) — ( dq ^\ + — = ^ q 

dr \Sdpr d ~ 1 dr J r d ~ l Sdpr^ 1 

Specializing to the case of a logotropic equation of state p = A In p, we obtain 

— ( A dq \ ° q = ^ 
1 } dr {SdP 2 ^- 1 dr) r*- 1 S d pr d - l<1 ' 

This equation is similar to Eq. (J75|) and they coincide at the point of marginal stability. There- 
fore, the stability threshold is the same in the two approaches and the above results concerning 
the form of the perturbation apply. Finally, if we consider the linear dynamical stability of a 
steady solution of the logotropic Euler-Poisson system ©-(J7J), we obtain an eigenvalue equation 
of the form ()86|) where A£ is replaced by A 2 [22] ■ Once again, the point of marginal dynamical 
stability coincides with the point where the energy functional ceases to be a minimum and 
becomes a saddle point. Therefore, the conditions of linear and nonlinear dynamical stability 
for the Euler-Poisson system coincide. 

4.8 Self-similar solutions of the polytropic and logotropic Smolu- 
chowski-Poisson system 

When no stable equilibrium state exists, the polytropic Smoluchowski-Poisson system admits 
self-similar solutions describing a gravitational collapse. These solutions have been given in 
|To] for a polytropic index n > 0. As we shall see, the situation is different, and richer, when 
n < (this case includes in particular the logotropes n = —1). We restrict our analysis to a 
space of dimension d > 2 (the dimension d = 2 is critical and requires a particular and intricate 
treatment; see [T7j for isothermal systems). 

If we introduce the function s(r, t) = M(r, t)/r d , which scales like the density p(r, t), we can 
rewrite the polytropic Smoluchowski-Poisson system (fTU|) - PHj) as a single differential equation 

m 

(87) * = 6 (r* + *) (*l + i±l*) + ( r * + ds 



dt \ dr J \dr 2 r dr J \ dr 

where we have defined A = if 7 and introduced an effective temperature O = A/S^ . We note 
that this equation passes to the limit for n — — 1. We look for self-similar solutions of the form 

s(r,t) = p (t)S 



ro(t) 

where p (t) is proportional to the central density. The scaling solution (|8Bj) only holds in the 
central region of the cluster. If we assume that all the terms scale the same in Eq. (J87j) . which 
needs not be the case (see below), we find that 

l/n+l 

(89) O^L „ p l 

'0 
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We can then define the scaling radius by 



(90) r 



e x 

l-l/n 



In that case, the relation between the typical central density and the typical core radius is 

2n 



(91) Po ~ r o °> with a 



n — 1 

Substituting these relations in Eq. (|57|) . we find that 



< 92 > * 



where we have set x = r/r . This implies that (1/ ' pfy(dpo/dt) is a constant that we arbitrarily 
set equal to a. This leads to 

(93) Po(t) = -(tcoii -t)~\ 

a 

so that the central density becomes infinite in a finite time t co u. The scaling equation now reads 

(94) aS + xS' = [xS' + dS) l/n fs" + ^jpS'J + {xS' + dS)S. 

For x — > +oo, we have the scaling S(x) ~ x _a so that the density behaves as p ~ r _Q for 
r — > +oo. At t = t co i; the density profile is proportional to l/r a . 

However, we could be in a situation in which the collapse is dominated by the gravitational 
drift like in the case of cold systems where 6 = 0. In that case, the dynamical equation (|87j) 
reduces to 

3s ( ds \ 

(95) _ =! ^_ + (fa j,. 

This equation with = has been solved in [TZJ . The collapse evolution is self-similar and the 
relation between the typical central density and the typical core radius is 

2d 

(96) po ~ r > Wlth a = 

We expect that this "cold regime" will prevail over the regime where the diffusion and the 
gravitational terms scale the same way if (p )e=o (po)e^o, i-e. if 

2d 2n 

(97) -r— > 



d + 2 n-l 

Indeed, one expects on the basis of free energy arguments that the most natural evolution is 
the one which leads to the most efficient collapse of the core. This is obtained by selecting the 
largest value of a between Eqs. (l9*Tj) and (|9*5jl. 

Let us consider the range of validity of this inequality: 

(i) The case n > has been considered in ^H]- It is shown that self-similar solutions exist 
only for n > n 3 = dj[d — 2) (otherwise the system tends to a complete polytrope). We note 
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Figure 6: Scaling exponent of the self-similar collapse as a function of the polytropic index (the 
curve is represented for d — 3). 

that inequality (}9Tj) is never satisfied in that case so that the scaling exponent is a = 2nj (n — 1) 
independant on the dimension. The results of are therefore unaltered. 

(ii) For n < 0, the inequality (|97|) is equivalent to n c = —d/2 < n < where n c = —d/2 is a 
new critical index which does not seem to have been introduced before. In that case, the collapse 
of the polytropic SP system is similar to the case = (cold system) and the scaling exponent 
is a = 2d/(d + 2). The corresponding scaling profile is known analytically in an implicit form 
[T7j . Finally, for n < n c = —d/2, the scaling exponent is a = 2n/(n — 1). In particular, for 
the logotropes with n = — 1 in d > 2, the scaling exponent is a = 2d/(d + 2). Therefore, the 
density profile of a collapsed logotrope scales as p ~ r -2d/(d+2) w ]jfl e ^he density profile of an 
equilibrium logotrope scales as p ~ r _1 for r — > +oo (see Sec. I4.4|) . This situation differs 
from the case of isothermal and polytropic systems with n > where the collapse exponent 
a = 2n/{n — 1) is the same as the exponent controlling the decay of the equilibrium density 
profile at large distances [T7| Ho]. 

The above results are summarized in Fig. |U1 In Fig. we show the data collapse obtained 
by solving the dynamical equation (J87|) with n = — 1 (logotrope) in dimensions d = 3 and 
d = 5. The scaling exponent agrees with its predicted theoretical value ()96|) corresponding to a 
pressureless collapse at O = 0. The scaling profile is also in good agreement with the analytical 
profile at © = obtained in in implicit form. In Fig. |H1 we show the time evolution 
of the density profile during the collapse obtained by solving the dynamical equation (jS7|) in 
dimension d = 3 for n = — 1 > n c = —3/2 and n = —4 < n c = —3/2. The exponents are in 
good agreement with the theoretical predictions (pjj) and (j9*Tj) respectively. 

5 Conclusion 

In this paper, we have studied a special class of nonlinear mean-field Fokker-Planck equations 
introduced in jH], corresponding to a logotropic equation of state. This can be viewed as 
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r/r (t) 

Figure 7: Data collapse for n = — 1 obtained by solving the dynamical equation (J%7|) . Upper 
curve: d = 3, a = 2d/(d + 2) = 6/5. Lower curve: d = 5, a = 2d/(d + 2) = 10/7; the curve has 
been shifted by an arbitrary factor. The invariant profile (dashed line) corresponds to = 0. 
It is known analytically in implicit form [T7] . 




Figure 8: The density profile (divided by Sd) is plotted for successive times so that the central 
density increases by a factor 2 at each time, starting from a central density 300. The dimension 
is d = 3. Upper curve: n = — 4 < n c = —d/2 so that a = 2n/(n — 1) = 8/5 (regular scaling); 
the figure has been shifted by a factor 100 for clarity. Lower curve: n = — 1 > n c so that 
a = 2d/(d + 2) = 6/5 (scaling at 6 = 0). The dashed lines are the theoretical slopes. 
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a limiting form of polytropic equation of state with 7 = 0. In the language of generalized 
thermodynamics, this corresponds to a Tsallis distribution with q = 0. Special attention has 
been given to the logotropic Smoluchowski-Poisson system where we have studied the steady 
states, their stability and, when unstable, the resulting gravitational collapse. In particular, 
we have shown that for polytropic indices — d/2 <n<0ind>2 (including the logotropes 
n = —1), the dynamics below the critical effective temperature G < @ c coincides with the 
dynamics at = 0. This is similar to the coarsening dynamics of spin systems strictly below 
the ferromagnetic critical temperature, which is controlled by the zero temperature fixed point 
|3*T] . This work completes previous investigations of the authors who studied the generalized 
Smoluchowski-Poisson system for an isothermal equation of state [HUE!, a polytropic equation 
of state with n > [Ej, a Fermi-Dirac equation of state |TH| and a Bose statistics in velocity 
space Our results can also have applications for the chemotactic aggregation of bacterial 
populations in biology that is governed by similar drift- diffusion equations [T%1 |2T| l2llj. 
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